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Abstract 

The slow drift (with speed e) of a parameter through a pitchfork bifurcation point, 
known as the dynamic pitchfork bifurcation, is characterized by a significant delay of 
the transition from the unstable to the stable state. We describe the effect of an addi- 
tive noise, of intensity a, by giving precise estimates on the behaviour of the individual 
paths. We show that until time after the bifurcation, the paths are concentrated in 
a region of size a/e^^^ around the bifurcating equilibrium. With high probability, they 
leave a neighbourhood of this equilibrium during a time interval [y/e, C\/e|logcr| ], after 
which they are likely to stay close to the corresponding deterministic solution. We 
derive exponentially small upper bounds for the probability of the sets of exceptional 
paths, with explicit values for the exponents. 
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1 Introduction 

Physical systems are often described by ordinary differential equations (DDEs) of the form 

n T 

^ = /(x,A), (1.1) 



where x is the state of the system, A a parameter, and s denotes time. The model (1.1) 
may however be too crude, since it neglects all kinds of perturbations acting on the system. 
We are interested here in the combined effect of two perturbations: a slow drift of the 
parameter, and an additive noise. 

A slowly drifting parameter \ = es, (with e <C 1), may model the deterministic change 
in time of some exterior influence, such as the climate acting on an ecosystem or a magnetic 
field acting on a ferromagnet. Obviously, nontrivial dynamics can only be expected when 
A is allowed to vary by an amount of order 1, and thus the system has to be considered 
on the time scale £~^. This is usually done by introducing the slow time t = es, which 
transforms ( |1.1D into the singularly perturbed equation 

n T 

e- = f{x,t). (1.2) 

It is known that solutions of this system tend to stay close to stable equilibrium branches 
of / |Gi] , Ti|, see Fig. |l^. New, and sometimes surprising phenomena occur when such an 
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Figure 1. Solutions of the slowly time-dependent equation (|l.2[ ) represented in the (i, x)- 
plane. (a) Stable case: A stable equilibrium branch x*{t) attracts nearby solutions xf°^. 
Two solutions with different initial conditions are shown. They converge exponentially fast 
to each other, as well as to a neighbourhood of order e of x*{t). (b) Pitchfork bifurcation: 
The stable equilibrium x = becomes unstable at t = (broken line) and expels two 
stable equihbrium branches ±x*{t). A solution xf°^ is shown, which is attracted hy x = 0, 
and stays close to the origin for a finite time after the bifurcation. This phenomenon is 
known as bifurcation delay. 



equilibrium branch undergoes a bifurcation. These phenomena are usually called dynamic 
bifurcations \ BenH H- In the case of the Hopf bifurcation, when the equilibrium gets unstable 



while expelling a stable periodic orbit, the bifurcation is substantially delayed: solutions 
of (|1.2D track the unstable equilibrium (for a non-vanishing time interval in the limit 



e ^ 0) before jumping to the limit cycle [ph] , Ne|. A similar phenomenon exists for the 
dynamic pitchfork bifurcation of an equilibrium without drift, the simplest example being 
f{x,t) = tx — (Fig- lllb)- The delay has been observed experimentally, for instance, in 



lasers [|ME[| and in a damped rotating pendulum [BK|. 

These phenomena have the advantage of providing a genuinely dynamic point of view 
for the concept of a bifurcation. Although one often says that a bifurcation diagram (rep- 
resenting the asymptotic states of the system as a function of the parameter) is obtained 
by varying the control parameter A, the impatient experimentalist taking this literally may 
have the surprise to discover unstable stationary states of the system (s)he investigates. 
The asymptotic state of the system (|1.1[) with slowly varying parameter A(es) = A(t) may 
depend not only on the initial condition (xq, to), but also on the history of variation of the 
parameter {X{t)}t^to- 



The perturbation of (1.1) by an additive noise can be modeled by a stochastic differ- 



ential equation (SDE) of the form 

dxs = f{xs,X)ds + adWs, (1.3) 

where Wg denotes the standard Wiener process, and a measures the noise intensity. A 
widespread approach is to analyse the probability density of Xg , which satisfies the Fokker- 
Planck equation. In particular, if — / can be written as the gradient of a potential function 
F, then there is a unique stationary density p{x, A) = e~^(^'^y'^ /N, where is the 
normalization. This formula shows that for small noise intensity, the stationary density is 
sharply peaked around stable equilibria of /. 

Unfortunately, the term "dynamical bifurcation" is used in a different sense in the context of random 
dynamical systems, namely to describe a bifurcation of the family of invariant measures as opposed to a 



"phenomenological bifurcation", see for instance |Ar 
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That method has, however, two major hmitations. The first one is that the Fokker- 
Planck equation is difficult to solve, except in the linear and in the gradient case. The 
second limitation is more serious: the density gives no information on correlations in 
time, and even when the density is strongly localized, individual paths can perform large 
excursions. This is why other approaches are important. A classical one is based on the 



computation of first exit times from the neighbourhood of stable equilibria [FW, FJ|. 



The effect of bifurcations has been studied more recently by methods based on the 



concept of random attractors [ CF94 , Schm , Ar |. In particular, Crauel and Flandoli showed 



that according to their definition, "Additive noise destroys a pitchfork bifurcation" |CF95]. 
The physical interpretation of random attractors is, however, not straightforward, and 
alternative characterizations of stochastic bifurcations are desirable. In the same way a 
slowly varying parameter helps our understanding of bifurcations in the deterministic case, 
it can provide a new point of view in the case of random dynamical systems. 

Let us consider the combined effect of a slowly drifting parameter and an additive noise 
on the ODE We will focus on the case of a pitchfork bifurcation, where the questions 

How does the additive noise affect the bifurcation delay? and Where does the path go after 
crossing the bifurcation point? are of major physical interest. The situation of the drift 



term / in (1.3) depending explicitly on time is considerably more difficult to solve than 
the autonomous case, and thus much less understood. One can expect, however, that a 
slow time dependence makes the problem accessible to perturbation theory, and that one 
may take advantage of techniques developed to study singularly perturbed equations such 



as (O). With X = £s, Equation (|1.3|) becomes 

dxs = f{xs,£s)ds + adWs. (1.4) 

If we introduce again the slow time t = es, the Brownian motion is rescaled, resulting in 
the SDE 

dxt = -fixt,t)dt + ^dWt. (1.5) 
£ ye 

Our analysis of ( |l.5D is restricted to one-dimensional x. The noise intensity a should be 
considered as a function of e. Indeed, since we now consider the equation on the time 
scale e~^, a constant noise intensity would lead to an infinite spreading of trajectories as 
e ^ 0. In the case of the pitchfork bifurcation, we will need to assume that a <^ ^fe. 

Various particular cases of equation ( |1.5D have been studied before, from a non-rigorous 
point of view. In the linear case f{x^ A) = Ax, the distribution of first exit times was 
investigated and compared with experiments in ]TM| , ^MC| , SIIA|, while pL[| derived a 



formula for the last crossing of zero. In the ..^{^=}^- x\ iJ^tudied the 



dependence of the delay on e and o numerically, while [Ku] considered the associated 
Fokker-Planck equation, the solution of which she approximated by a Gaussian Ansatz. 

In the present work, we analyse ( |1.5| ) for a general class of odd functions /(x. A) under- 
going a pitchfork bifurcation. We use a different approach, based on a precise control of 
the whole paths {xs}to^s^t of the process. The results thus contain much more information 
than the probability density. It also turns out that the technique we use allows to deal 
with nonlinearities in quite a natural way. Our results can be summarized in the following 
way (see Fig. ^): 

• Solutions of the deterministic equation ( |1.2D starting near a stable equilibrium branch 
of / are known to reach a neighbourhood of order e of that branch in a time of order 
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Figure 2. A typical path xt of the stochastic differential equation (1.5) near a pitchfork 
bifurcation. We prove that with probability exponentially close to 1, the path has the 
following behaviour. For ^ t ^ y^, it stays in a strip B{h) constructed around the 
deterministic solution with the same initial condition. After t = -^/e, it leaves the domain 
2? at a random time t — tx>, which is typically of the order y^ejlogfTl. Then it stays 
(up to times of order 1 at least) in a strip A^{h) constructed around the deterministic 
solution a;f''''^ starting at time t on the boundary of T>. The widths of B{h) and A^{h) 
are proportional to a parameter h satisfying a ^ h ^ y/e. 



e|loge|. We show that the paths of the SDE ( |1.5D with the same initial condition 
are typically concentrated in a neighbourhood of order a of the deterministic solution 
(Theorem 2.3). 



• A particular solution of the deterministic equation (1^) is known to exist in a neigh- 
bourhood of order e of each unstable equilibrium branch of /. Paths that start in a 
neighbourhood of order a of this solution are likely to leave that neighbourhood in a 
time of order e|loge| (Theorem 2.5). 

• When a pitchfork bifurcation occurs at x = 0, t = 0, the typical paths are concentrated 
in a neighbourhood of order a je^l'^ of the deterministic solution with the same initial 
condition up to time \/e (Theorem ^.S]) . 

• After the bifurcation point, the paths are likely to leave a neighbourhood of order \/t 
of the unstable equilibrium before a time cy^e|logcj| (Theorem p^ ). 

• Once they have left this neighbourhood, the paths remain with high probability in a 
region of size o jyft around the corresponding deterministic solution, which approaches 
a stable equilibrium branch of / like eji?!'^ (Theorem 2.10 ). 

These results show that the bifurcation delay, which is observed in the dynamical sys- 
tem ([1.2|), is destroyed by additive noise as soon as the noise is not exponentially small. 
Do they mean that the dynamic bifurcation itself is destroyed by additive noise? This is 
mainly a matter of definition. On one hand, we will see that independently of the initial 
condition, the probability of reaching the upper, rather than the lower branch emerging 
from the bifurcation point, is close to \. The asymptotic state is thus selected by the 
noise, and not by the initial condition. Hence, the bifurcation is destroyed in the sense of 
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| CF98 |. On the other hand, individual paths are concentrated near the stable equilibrium 
branches of /, which means that the bifurcation diagram will be made visible by the noise, 
much more so than in the deterministic case. So we do observe a qualitative change in 
behaviour when A changes its sign, which can be considered as a bifurcation. 

The precise statements and a discussion of their consequences are given in Section 
In Section we analyse the motion near equilibrium branches away from bifurcation 



points. The actual pitchfork bifurcation is discussed in Section 2.3. A few consequences 



are derived in Section |2.4| . Section y contains the proofs of the first two theorems on the 
motion near nonbifurcating equilibria, while the proofs of the last three theorems on the 
pitchfork bifurcation are given in Section ^. 
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2 Statement of results 
2.1 Preliminaries 

We consider nonlinear Ito SDEs of the form 

dxt = -f{xt,t)dt + -^dWt, xto=xo, (2.1) 
e ye 

where {Wt}t^to is the standard Wiener process on some probability space {Q,T,F). Ini- 
tial conditions xq are always assumed to be square-integrable with respect to P and in- 
dependent of {Wt}t^to- All stochastic integrals are considered as Ito integrals, but note 
that Ito and Stratonovich integrals agree for integrands depending only on time and iv. 
Without further mentioning we always assume that / satisfies the usual (local) Lipschitz 
and bounded-growth conditions which guarantee existence and (pathwise) uniqueness of a 
(strong) solution {xt}t of (|2.l| ). Under these conditions, there exists a continuous version 
of {xt}t- Therefore we may assume that the paths uj i— > xt{uj) are continuous for P-almost 
all UJ £ Q. 

We introduce the notation W^o,xo fQj- gf ^j^g process {xt}t:^to^ starting in xq 

at time to, and use E*"'^" to denote expectations with respect to P*0'^o. Note that the 
stochastic process {xt}t^to is an (inhomogeneous) Markov process. We are interested in 
first exit times of xt from space-time sets. Let ^ C M x [to,^!] be Borel-measurable. 
Assuming that A contains (xo,to)) we define the first exit time of {xt,t) from A by 

TA = mi{t€[to,ti]: {xt,t)^A}, (2.2) 

and agree to set t_4^{u>) = oo for those uj £ fl which satisfy (xj(w), t) G A for all t G [to, h]. 
For convenience, we shall call the first exit time of xt from A. Typically, we will 
consider sets of the form A = {{x,t) G M x [to,ti]: 51 (i) < x < g2{t)} with continuous 
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functions gi < 52- Note that in this case, r_4 is a stopping timeQ with respect to the 
canonical filtration of {Q,J^,F) generated by {xt}t^to- 

Before turning to the precise statements of our results, let us introduce some notations. 
We shall use 

• [y] for y ^ to denote the smallest integer which is greater than or equal to y, and 

• yy z and y Az to denote the maximum or minimum, respectively, of two real numbers 
y and z. 

• By g{u) = 0{u) we indicate that there exist 6 > and K > such that g{u) ^ Ku for 
all u G [0, 6], where 6 and K of course do not depend on e or a. Similarly, g{u) = o{l) 
is to be understood as liniu^o g{u) = 0. From time to time, we write g{u) = Ot{1) to 
indicate that choosing a priori a sufficiently small T allows to make the corresponding 
term arbitrarily small for all u from some T-dependent interval. 

Finally, let us point out that most estimates hold for small enough e only, and often only 
for P-almost all uj £ 0,. We will stress these facts only when confusion might arise. 



2.2 Nonbifurcating equilibria 

We start by considering the nonlinear SDE ( |2.1| ) in the case of / admitting a nonbifurcating 
equilibrium branch. We will assume that there exists an interval / = [0, T] or [0, 00) such 
that the following properties hold: 

• there exists a function x* : / ^ M , called equilibrium curve, such that 

/(x*(t),t) = VtGl; (2.3) 

• / is twice continuously differentiable with respect to x and t, with uniformly bounded 
derivatives, for all t G I and all x in a neighbourhood of x*{t); 

• the linearization of / at x*{t), defined as 

a{t) = dj{x^{t),t), (2.4) 
is bounded away from zero, that is, there exists a constant oq > such that 

\a{t)\ ^ao Vt G /. (2.5) 
In the deterministic case a = 0, the following result is known (see Fig. |l|a): 



Theorem 2.1 (Deterministic case [ Fi , |Gij| ) . Consider the equation 



e^ = f{xt,t). (2.6) 

There are constants £0, cq, ci > 0, depending only on f , such that for < e ^ eo; 

• ( |2.6| ) admits a particular solution x^'^^ such that 

-x''(t)K cie yt£l; (2.7) 

• "if \xo — x*{0)\ ^ Co and a{t) ^ —a^ for all t £ I (that is, when x* is a stable equilib- 
rium), then the solution xf^^ of ( p.6D with initial condition Xq'^* = xq satisfies 

l^det _ jdot| ^ _ jdet| ^-aot/2e ^ j (2.8) 



■^For a general Borel-measurable set A, the first exit time ta is still a stopping time with respect to the 
canonical filtration, completed by the null sets. 
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Remark 2.2. The particular solution x'^'^^ is often called a slow solution or adiabatic 
solution of equation ( |2.6| ). It is not unique in general, as suggested by (^^). 



We return now to the SDE (2.1) with a > 0. We need no additional assumption on 
a in this section. However, the results are only interesting when a = Os{l)- Let us first 
consider the stable case, that is, we assume that a{t) ^ —ao < for all t £ I. We assume 
that at t = 0, xt starts at some (deterministic) xq sufficiently close to x*{0). Theorem 

xf"^^ with the same initial condition Xq^^ = xq 



2A\ tells us that the deterministic solution x?*^* with the same initial condition ""'^^^ 



reaches a neighbourhood of order e of x*(t) exponentially fast. 

We are interested in the stochastic process yt = Xf — xf^^, which describes the deviation 
due to noise from the deterministic solution x'^*^*. It obeys the SDE 

dyt = - [f{xf'' + yt, t) - f{xf\t)] dt + ^ dWt, yo = 0. (2.9) 
£ ye 

We will prove that yt remains in a neighbourhood of with high probability. It is instruc- 
tive to consider first the linearization of ( ^.9| ) around y = 0, which has the form 

dy° = -a(t)yO dt + ^ dWt, (2.10) 
e ye 

where 

a{t) = dj{xf'\ t) = a{t) + 0{e) + O{\xo - x*(0)| e-"«*/2^). (2.11) 

Taking e and \xo — x*(0)| sufficiently small, we may assume the existence of constants 
a+ ^ a_ > such that — a+ ^ a{t) ^ — a_ for all t £ I. The solution of ( |2.1C1| ) with 
arbitrary initial condition y^ is given by 

yO = yO ^m/e +^ /"* ^^^^^ ^(^^ ^ /■* (2.12) 

Jo Js 

where we write a{t, 0) = a{t) for brevity. Note that a{t, s) ^ — a_(t — s) whenever t ^ s. 
If j/q has variance fo ^ 0, then y^ has variance 

v{t) = vo e2°(*)/^ +— /"* e2°(*'^)/^ ds. (2.13) 
^ Jo 

Since the first term decreases exponentially fast, the initial variance vq is "forgotten" as 
soon as e^"^*)/'^ is small enough, which happens already for t > 0{e\loge\). For t/q = 0, 
(2.12) implies in particular that for any 5 > 0, 

P°'°{l2/?l ^5} ^e-'^'/2-W, (2.14) 

and thus the probability of finding y^, at any given t G I, outside a strip of width much 
larger than y^2v{t) is very small. 

Our first main result states that the whole path {x^jo^s^t of the solution of the non- 
linear equation ( |2.1| ) lies in a similar strip with high probability. We only need to make 
one concession: the width of the strip has to be bounded away from zero. Therefore, we 
define the strip as 

Bs{h) = {{x,t) eR X I:\x- < hy^}, (2.15) 
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where 



2|a(0)| 



+ - 



(2.16) 



(T^C can be interpreted as the variance (2.13) of the process ( |2.12 ) starting with initial 
variance vq = (T^/(2|a(0)|). We shall show in Lemma that 



2\a{t)\ 



+ Oie) + 0(\xo-x*{0)\e 



-aot/2e 



(2.17) 



Let Tg^(;i) denote the first exit time of Xf from Bs{h). 



Theorem 2.3 (Stable case). There exist eq, do and ho, depending only on f, such that 
for < e ^ £o, h ^ ho and \xo — x*(0)| ^ do, 



{rB4h) <t}^ C{t, e) exp{-i^ [1 - 0{e) - 0{h)\ }, (2.18) 



where 



Cit,e) 



\ait)\ 



+ 2. 



(2.19) 



The proof, given in Section is divided into two main steps. First, we show that 
an estimate of the form ( |2.18| ), but without the term 0{h), holds for the solution of the 
linear equation ( p. 101 ). Then we show that whenever \y^\ < (^[s) for ^ s ^ i, one 
almost surely also has \ys\ < h{l + 0{h))y/(^[s) for ^ s ^ t. 

Remark 2.4. The result of the preceding theorem remains true when l/2|a(0)| in the 
definition ( 2.16| ) of C(t) is replaced be an arbitrary (q, provided Co > 0. The terms 0{-) may 
then depend on (q. Note that ({t) and a'^v{t) are both solutions of the same differential 
equation sz' = 2a(t)z + 1, with possibly different initial conditions. If xq — x*(0) = 0(e), 
({t) is an adiabatic solution (in the sense of Theorem p.l| ) of the differential equation, 
staying close to the equilibrium branch z* = l/|2a(t)|. 

The estimate ( p. 18 ) has been designed for situations where o" <C 1, and is useful for 
(J h 1. We expect the exponent to be optimal in this case, but did not attempt to 
optimize the prefactor C(t,e), which leads to subexponential corrections. If we assume, 
for instance, that cr = e'^, g > 0, and take h = with Q < p < q, ( 2.18 ) can be written as 

F°'^n^B.(h) <t]^ (t + e2)exp{-^-^[l - 0{e) - Oie^ - 0(e2(«-^)|log£|)] }. 

(2.20) 

The t-dependence of the prefactor is to be expected. It is due to the fact that as time 
increases, the probability of xi escaping from a neighbourhood of x^^^ also increases, but 
very slowly if a is small. The estimate ( p. 18 ) shows that for a fraction 7 of trajectories to 
leave the strip Bs{h), we have to wait at least for a time t^ given by 



which is compatible with results on the autonomous case. 



(2.211 
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Let us now consider the unstable case, that is, we now assume that the linearization 



a{t) = dxf{x*{t),t) satisfies a{t) ^ ao > for all tel. Theorem 2A shows the existence of 
a particular solution xf"^^ of the deterministic equation (^) such that \xf'^^ — x*{t)\ ^ cie 
for all tel. We define a{t) = dxf{xf^\t) = a{t) + 0{e) > and a{t) = a{s) ds. 

The linearization of around xf^^ again admits a solution of the form ( ^.12| ) . In 
this case, however, the variance ( p. 13 ) grows exponentially fast, and thus one expects the 



probability of xt remaining close to xf'^^ to be small. This is the contents of the second 
main result of this section. We introduce the set 

B^{h) = |(x,t) gM X/: |x-xf*| < \ \ (2.22) 
I ^/2a{t)] 

and the first exit time Tj^^i^h) of xt from B^i{h). 

Theorem 2.5 (Unstable case). There exist eq and ho, depending only on f , such that 
for all h ^ a A Kq, all e ^ Sq and all xq satisfying (xqjO) G Bu{h), we have 

where k = f^{l - 0{h) - 0{e)) . 



The proof, given in Section 3.2, is based on a partition of the interval [0,t] into small 
intervals, and a comparison of the nonlinear equation with its linearization on each interval. 

This result shows that xt is unlikely to remain in B^^h) as soon as t > ea'^/h^. A 
major limitation of ( |2.23| ) is that it requires h ^ a. Obtaining an estimate for larger h is 
possible, but requires considerably more work. We will provide such an estimate in the 



more difficult, but also more interesting case of the pitchfork bifurcation, see Theorem 2.£ 
below. 



2.3 Pitchfork bifurcation 

We now consider the SDE ( p.l| ) in the case of / undergoing a pitchfork bifurcation. We 
will assume that 

• / is three times continuously differentiable with respect to x and t in a neighbourhood 
^fo of (0,0); 

• f{x,t) = -fi-x,t) for all (x,t) G Mq] 

• / exhibits a supercritical pitchfork bifurcation at the origin, i.e. 

9,/(0,0) = 0, at,/(0,0)>0 and 9,,,/(0, 0) < 0. (2.24) 

The assumption that / be odd is not necessary for the existence of a pitchfork bifur- 
cation. However, the deterministic system behaves very differently if x = is not always 
an equilibrium. The most natural situation in which f{0,t) = for all t is the one where 
/ is odd. 

By rescaling x and t, we may arrange that 5jx/(0,0) = 1 and dxxxfiO,0) = —6 as in 
the standard case f{x,t) = tx — x^. This implies in particular that the linearization of / 
at X = satisfies 

a{t) = dxf{Q,t) = t + 0{t'). (2.25) 
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A standard result of bifurcation theory | GH , |j| states that under these assumptions, there 
is a neighbourhood M C A/q of (0, 0) in which the only solutions of /(x, t) = are the line 
X = and the curves 

x = ±x*{t), x*{t) = \ft[l + Ot{l)\, t^O. (2.26) 

If M is small enough, the equilibrium x = is stable for t < and unstable for t > 0, 
while X = zizx*{t) are stable equilibria with linearization 

a*{t) = dJix*{t),t) = -2t[l + Ot{l)]. (2.27) 

The only solutions of dxf{x,t) = in M are the curves 

x = ±xit), x{t) = y/t/3[l + Otil)], t^O. (2.28) 

If / is four times continuously differentiable, the terms 0^(1) in the last three equations 
can be replaced by 0{t). 

We briefly state what is known for the deterministic equation 

e^ = fixt,t), (2.29) 

where we take an initial condition (xo,to) £ with xq > and to < 0, see Fig. |l|b. 
Observe that a{t, to) = /j* a{s) ds is decreasing for to < ^ < and increasing for t > 0. 

Definition 2.6. The bifurcation delay is defined as 

n(to) = inf{t > 0: a(t,to) > O}, (2.30) 
with the convention n(to) = oo if a{t,t()) < for all t > 0, for which a(t,to) is defined. 

One easily shows that n(to) is differentiable for to sufficiently close to 0, and satisfies 
limt(,_o- n(to) = and limt^^o- n'(to) = -1. 

Theorem 2.7 (Deterministic case). Let xf^'^ be the solution of ( ^.291) with initial con- 
dition xf^^ = xq- Then there exist constants Eq, cq, c\ depending only on f, and times 

ti = to + 0{e\logs\) 

t2 = n(ti) = n(to)-0(e|loge|) (2.31) 
t3 = n(to) + 0(e|loge|) 

such that, i/0 < xo ^ Co, < e ^ eo O'^d {xf'^^,t) G Af, 

fo < xf * ^ ciee'^^*'*!)/^ for ti ^ t ^ t2 ,^ 
||x^^* - x*(t)| ^ cie forf^ts. 

The proof is a straightforward consequence of differential inequalities, see for instance 



| Ber , Propositions 4.6 and 4.8]. 

We now consider the SDE ( |2.1| ) for o" > 0. The results in this section are only inter- 
esting for a = o{y/e), while one of them (Theorem |2.9| ) requires a condition of the form 
cr|log(Tp/^ = 0{^/e) (where we have not tried to optimize the exponent 3/2). 
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Let us fix an initial condition {xtg,to) £ M witli to < 0. For any T G (0, ItoD) we 



can apply Theorem 2.3 on the interval [to, —T] to show that \x-t\ is likely to be of order 
a^^^ + cie e°'(~'^'^'^y^ for any S > 0. We can also apply the theorem for t > T to show that 
the curves zizx*(t) attract nearby trajectories. Hence there is no limitation in considering 
the SDE ( |2.lD in a domain of the form |x| ^ d, |t| ^ T where d and T can be taken 
small (independently of e and a of course!), with an initial condition x-t = xq satisfying 
|xo| ^ d. 

We first show that xt is likely to remain small for — T ^ t ^ y/e. Actually, it turns 
out to be convenient to show that xt remains close to the solution xo e"^*'""^)/'^ of the 
linearization of ( |2.29| ). We define the "variance- like" function 

C(t) = — ^ e2°(*'-^)/- +- f e2"(*'^)/- ds. (2.33) 
2|o(-r)| e J^T 



We shall show in Lemma |4.2| that for sufficiently small e, there exist constants c± such 
that 

Y\ ^ for -T ^ t ^ (2.34) 

^ C{t) ^ for -y/e ^t^^/e. (2.35) 

The function C(t) is used to define the strip 

B{h) = {{x,t) e [-d,d] X [~T,^/I]■. |x-xoe"(*'-^)/=| < h^/C(t)}. (2.36) 
Let Tg(/i) denote the first exit time of xt from B{h). 

Theorem 2.8 (Behaviour for t ^ \/e)" There exist constants andhQ, depending only 
on f , T and d, such that for < e ^ Eq, h ^ ho^/e, |xo| ^ h/e^^^ and —T ^ t ^ ^/e, 

IP~^'"''{TB(h)<t}^C(t,e)exp|-i^ l-r{e)-o(^^'^ | (2.37) 

where 

Cft,e)= l°"--l + °<^' . (2.38) 

and with r(e) = ©(e) for -T ^ t ^ -^/e, and r(e) = ©(V^) for -^/e ^ti^^/e. 



The proof (given in Section O) and the interpretation of this result are very close in 
spirit to those of Theorem |2.3| . The only difference lies in the kind of e-dependence of the 
error terms. The estimate ( |2.37] ) is useful when a <^ h <^ y^, and shows that the typical 
spreading of paths around the deterministic solution will slowly grow until t = ^/e, where 
it is of order a/e^^^, see Fig. |2[ 

Let us now examine what happens for t ^ ^/e. We first show that xt is likely to leave 
quite soon a suitably defined region D containing the line x = 0. The boundary of T> is 
defined through a function x(t), which can be chosen somewhat arbitrarily, but should lie 
between x(t) and x*(t), in order to simplify the analysis of the dynamics after xt has left 
T>. A convenient definition is 

x{t) = VXx*{t), (2.39) 
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where A is a free parameter. We need to assume, however, that A G We now define 

V = {{x,t) G [-d,d] X [^/e,T]■. \x\ < (2.40) 
Note that V has the property that for all (x,t) G V with a; / 0, 

-f(x,t) ^ Ka(t) with k; = I-A-cit(I). (2.41) 

X 

Let Tx) denote the first exit time of xt from D. 

Theorem 2.9 (Escape from V). Let (xo,to) S V and assume thata\logaf/'^ = 0{^/E). 
Then for to ^t^T, 

p-.{.. > *} < c„iWv^M^(i + ^) ^ „ , (2.42) 

a \ £ J _ Q-2Ka{t,to)/e 

where Cq > is a (numerical) constant. 

The proof of this result (given in Section is by far the most involved of the present 
work. We start by estimating, in a similar way as in Theorem 2^, the first exit time 
from a strip S of width slightly larger than u/y^ a{s). The probability of returning to zero 
after leaving S can be estimated; it is small but not exponentially small. However, the 
probability of neither leaving P nor returning to zero is exponentially small. This fact can 
be used to devise an iterative scheme that leads to the exponential estimate ( p.42| ) . 

We point out that for every subset V C V, we have P*0'^o{ri,/ ^ t} ^ P*o.^o{^^ ^ 
and thus ( p. 42 ) still provides an upper bound for the first exit time from smaller sets. 

Let us finally consider what happens after the path has left P at time t = tx>. One 
can deduce from the definition ( 2.39| ) of x{t) that for ^/e ^ t and |x| ^ x{t), 

d,,f{x,t) f^d{t) = d:^f{x{t),t) f^-rjait) with r/ = 3A - 1 - Ot(1). (2.43) 

Let x^''*'^ denote the solution of the deterministic equation ( p. 29)) starting in x{t) at 
time T (the case where one starts at —x{t) is obtained by symmetry). We shall show in 
Proposition 4.11 that x^^ always remains between x(t) and x*{t), and approaches x*{t) 
according to 



= x*{t) - 0[^) - 0{V^e-^^(''^y^). (2.44) 
Moreover, deterministic solutions starting at different times approach each other like 

^ xf - ^ (x^"*'v^ - x(r)) e-''"^*'")/^ Vt G [r, T]. (2.45) 



The linearization of / at xf^^''^ satisfies 



a^t) = a,/(xf*'^^) = a*(t) + +0(te-''°(*'^)/^). (2.46) 

For given r, we construct a strip A'^{h) around x*^^*'"^ of the form 

A^ih) = {{x,t): T ^ t ^r,|x-xf*'^| < h^/C{t)}, (2.47) 
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where the function C^(^) is defined by 



at) 



1 



2|a(r) 



e2°^(*'^)/^ ds, a-{t, s) = I' a^u) 



dn, 



and satisfies 



at) 



1 



2|a*(t)| 



t3 



(2.48) 



(2.49) 



cf. Lemma 4.12| . Let t^^^/j) denote the first exit time of Xt from A'^{h). 

Theorem 2.10 (Approach to x*). There exist constants £q and ho, depending only on 
f, T and d, such that for < e ^ eo? ^ < h^r and t ^ t ^ T , 



\)T,x{t) 



{TA^k) <t}< C^it,E) exp{-^^ [l - (Die) - } 



where 



a*(s) ds 



+ 2. 



(2.50) 



(2.51) 



The proof is given in Section |4.4| . This result is useful for a <^ h <^ t, and shows that 
the typical spreading of paths around x^^^'^ is of order a j^fi^ see Fig. |2|. 



2.4 Discussion 

Let us now examine some of the consequences of these results. First of all, they allow to 
characterize the influence of additive noise on the bifurcation delay. In the deterministic 
case, this delay is defined as the first exit time from a strip of width e around x = 0, 
see Theorem 2.7. A possible definition of the delay in the stochastic case is thus the first 
exit time t^^^'^^ from a similar strip. An appropriate choice for the width of the strip is 
x(v^) = C(e^/^), since such a strip will contain B(h) for every admissible /i, and the part 
of the strip with t ^ yfe will be contained in T). Theorems 2^ and 2^ then imply that if 
t ^ a/F, 



-T,xo|^delay < ^} ^ C( V^, e) e" 



.0(£/<72) 



-T,xo|^delay ^ ^| ^ C^i^t)^^ 



|logo-| 



1 + 



a 



-2Ka(t,v^)/e 



(2.52) 
(2.53) 



If we choose t in such a way that a(t, \fe) = ce\loga\ for some c > 0, the last expression 
reduces to 

p-T,xo|^delay ^ ^| ^ O ((7^'^-^ |log a] 2) , (2.54) 

which becomes small as soon as c > The bifurcation delay will thus lie with over- 

whelming probability in the interval 

[V^,0{^/e\loga\)]. (2.55) 



Theorem 2.10| implies that for times larger than 0{y^ e|log a\ ), the paths are unlikely to 
return to zero in a time of order 1. The wildest behaviour of the paths is to be expected 
in the interval ( 2.55| ), because a region of instability is crossed, where dxf > 0. 

Our results on the pitchfork bifurcation require a ^ ^/£, while the estimate ( p. 55] ) 
is useful as long as a is not exponentially small. We can thus distinguish three regimes, 
depending on the noise intensity: 
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a ^ y/e: A modification of Theorem 2^ shows that for t < —a, the typical spreading 
of paths is of order a/^^/\t\. Near the bifurcation point, the process is dominated 
by noise, because the drift term / ~ — is too weak to counteract the diffusion. 
Depending on the global structure of /, an appreciable fraction of the paths might 
escape quite early from a neighbourhood of the bifurcation point. In that situation, 
the notion of bifurcation delay becomes meaningless. 

g-i/eP ^ ^ for some p < 1: The bifurcation delay lies in the interval ( 2.55| ) with 



high probability, where Y^e|log(T| ^ still "microscopic 

^ e-^l^ for some K > The noise is so small that the paths remain concentrated 
around the deterministic solution for a time interval of order 1. The typical spreading 
is of order a\J C{t), which behaves like cre"^*^/^ /e^^'^ for t ^ y/e, see Lemma [4.2| . Thus 
the paths remain close to the origin until a{t) ~ e|log(j| ^ K. If e|log(j| > a(n(io)) = 
|Q;(io)|, they follow the deterministic solution which makes a quick transition to x*{t) 

at t = n(to). 

The expression ( |2.55 ) characterizing the delay is in accordance with experimental results 



| TM , SMC I, and with the approximate calculation of the last crossing of zero [|JL|] . The 



numerical results in |Ga|, which are fitted, at e = 0.01, to r'^^^^y ~ ^0-105 weak noise 



and r'^^^'^y ~ e ^^^'^ for strong noise, seem rather mysterious. Finally, the results in [Ku], 



who approximates the probability density by a Gaussian centered at the deterministic 
solution, can obviously only apply to the regime of exponentially small noise. 

Another interesting question is how fast the paths concentrate near the equilibrium 
branches zizx*{t). The deterministic solutions, starting at x(to) at some time to > 0, all 
track x*{t) at a distance which is asymptotically of order e/t^/^. Therefore, we can choose 
one of them, say and measure the distance of xt from that deterministic solution. 

We restrict our attention to those paths which are still in a neighbourhood of the origin 
at time ^/e, as most paths are. We want to show that for suitably chosen ti E {^/£, t) and 
A G (0, t), most paths will leave V until time ti and reach a ^-neighbourhood of x^*^*'^ 
at time + A. Let us estimate 



TX) <ti, sup 

se[ri3+A,t] 



dot.v^ 



<(5j j (2.56) 

I Ue[rT,+A,tl J J 



The first term decreases roughly like cr ^ e '^"{*i'V^)/'^ and becomes small as soon as 
a{tl,^/e) » e|logcr|. The second summand is bounded above by 

const E^'^'^{l|,^<i^} ^^P{- J - 0(V^e-''"(-^+^'-^)/^)'] }}. (2.57) 

Therefore, 6 should be large compared to a/t and we also need that A is at least of order 
0(ye|logcj|). Then we see that after a time of order 0{y/£\loga\), the typical paths 
will have left D and, after another time of the same order, will reach a neighbourhood of 
^det,^^ which scales with a/t. 

Finally, we can also estimate the probability of reaching the positive rather than the 
negative branch. Consider Xg, starting in xq at time to < 0, and let t > 0. Without loss 
of generality, we may assume that xq > 0. The symmetry of / implies 

P*0'^o{xt ^ 0} = 1 - lp*o-^o|3 5 g [ip^^) . ^ 0}, (2.58) 
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and therefore it is sufficient to estimate the probabihty for Xg to reach zero before time 
zero, for instance. We Unearize the SDE (2.1) and use the fact that the solution x^ of the 
Unearized equation 



dx!^ 



-a(s)a;° ds + ^ dWs 
£ ye 



% = ^0 



(2.59) 



satisfies Xg ^ x^ as long as Xg does not reach zero. For the Gaussian process x^ we know 



P*»'^o{3s E [to,t) : x° = 0} = 2(l-P*«'^o{2;° > O}) = 1- ^ / 



(2.60) 



where u{t) = xq e°(*'*(')/^ /y^v{t,to) and v{t,to) denotes the variance of x^. For t = 0, 
u{0) is of order xo£^^^a~^ e~^°'^^^ ^o/^ ^ see Lemma 4^. Thus the probability in (|2.60|) is 
exponentially close to one for small e, and we conclude that the probability for xt to reach 
the positive branch rather than the negative one is exponentially close to 1/2. 



3 The motion near nonbifurcating equilibria 

In this section we consider the nonlinear SDE 



dxt = -f{xt,t)dt + ^dWt (3.1) 
£ ye 

under the assumptions 

• t E / = [0,r] or [0,oo); 

• there exists an equilibrium curve x* : I such that 

f{x*{t),t)=0 VtE/; (3.2) 

• there is a constant d > such that / is twice continuously differentiable with respect 
to X and t for \x — x*{t)\ ^ d and t £ I, with \dxxf{x, t)\ uniformly bounded by 2M > 
in that domain; 

• there is a constant gq > such that a{t) = dxf{x*{t),t) satisfies 

\a{t)\ ^ ao Vt E /. (3.3) 

We do not need any assumptions on u > 0, but our results are of interest only for a = ©£(1). 

In Section ^]l] we consider the stable case, corresponding to a{t) ^ — ao < for all 
t £ I. We first analyse the linearization of (|3.1D around a given deterministic solution. 



Proposition 3.3 shows that the solutions of the linearized equation are likely to remain 



in a strip of width hyQ{f) around the deterministic solution. Here C,{t) is related to 



the variance and will be analyzed in Lemma |3.l| . Proposition |3.6| allows to compare the 



trajectories of the linear and the nonlinear equation, and thus completes the proof of 



Theorem 2.3. 



In Section ^y^, we consider the unstable case, i. e. a{t) ^ ao > for all t £ I. 
Theorem |2.5| is equivalent to Proposition ^j^, which is again based on a comparison of 
solutions of the nonlinear equation (|3.1| ) and its linearization around a given deterministic 
solution. 
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3.1 Stable case 

We first consider tlie case of a stable equilibrium, that is, we assume that a(t) ^ —uq for 
ah t & I. We win assume that the stochastic process xt, given by the SDE (|3.lD, starts 



at time t = in xq- By Theorem 2.1, there exists a cq > such that the deterministic 
solution x'^'^^ of ( |2.6D with initial condition Xq'^* = xq satisfies 

|xf * - ^ 2cie + |xo - x*(0)| e-''«*/2e Vt e /, (3.4) 

provided |xo — x*(0)| ^ cq. We are interested in the stochastic process yt = xt — xf'^^, 
which describes the deviation due to noise from the deterministic solution x'^^^. It obeys 
an SDE of the form 

dyt = -[a{t)yt + b{yt,t)]dt + ^dWt, yo = 0, (3.5) 
£ ye 

where we have introduced the notations 

„det 



a{t) = ae{t)=d,f{xr\t) 
b{y, t) = b,{y, t) = /(xf * +y,t)- f{xf'\t) - a{t)y. 



(3.6) 



Taking e and |xo — x*(0)| sufficiently small, we may assume that there exists a constant 
J > such that Ix^'^* + y — x*{t)\ ^ d whenever \y\ ^ d. It follows from Taylor's formula 
that for all (y, t) G [—d, d] x /, 

\b{y,t)\^My^ (3.7) 
\a{t) - a{t)\ < M{2cie + |xo - x*(0)| e-""*/^^) (3_8) 

By again taking e and |xo — 2;*(0)| sufficiently small, we may further assume that there 
are constants a+ ^ a_ > ao/4 such that 

-a+ ^ a{t) ^ -a_ G /. (3.9) 

Finally, the relation a'{t) = dxtf{xf'^^,t) + dxxf{xf'^^,t)^f{xf^^,t) implies the existence of 
a constant C2 > such that 

\d'{t)\ ^ C2(l + |xo-x*(0)| . (3.10) 



Our analysis will be based on a comparison between solutions of ( |3.5| ) and those of the 
linearized equation 

dyO = -a(t)y° dt + ^ dWt, = 0. (3.11) 

Its solution is given by 

yt=^ / e°(*'^)/^dWs, a(t,s)= / a(^/)du. (3.12) 

Jo Js 

We will write 5(t, 0) = 5(t) for brevity. The Gaussian random variable y^ has mean zero 
and variance 

v{t) = — re2ffi(t,s)/e^^_ ^3_;^3) 
£ Jo 



16 



Note that ( |3.9D imphes that a{t,s) ^ — a_(t — s) whenever t ^ s, which imphes in 
particular, that v{t) is not larger than (T^/2a_. We can, however, derive a more precise 
bound, which is useful when e and e~°-ot/'i£ ^j-g gniall. To do so, we introduce the function 



at) 



2|a(0)| 



ds, where a{t) = a{t,0). 



(3.14) 



Note that v{t) ^ a'^(^[t), and that both functions differ by a term which becomes negligible 
as soon as t > ©(ellogej). The behaviour of C{t) is characterized in the following lemma. 



Lemma 3.1. The function ((t) satisfies the following relations for all t £ I. 

1 



2\a{t)\ 



+ O{e) + O{\xo- 



-aot/2e) 



2a 



2a- 



C'it) ^ ^ 



Proof: By integration by parts, we obtain that 

1 



Using IK 



and ( |3l0D we get 



-2a{t) 



* a'is) 



a(s)2 



,2a{t,s)/e 



ds. 



(3.15) 
(3.16) 
(3.17) 

(3.18) 



' -a' is) , 





,2a(t,s)/e 



2a_ {t-s)/e + \xo-X*{0)\ ^[_2s_ (t_,)_„„,/2]/^ ^ 



C2 C2 Xo 

^ — -e + 



x*m 



-aot/2e 



2al c?_ 2a- -ao/2'^ 
which proves ( 3.15| ). We now observe that ({t) is a solution of the linear ODE 

dC 



dt 



:(2a(t)C + l), 



C(0) 



2|a(0)| ' 



(3.19) 



(3.20) 



Since C{t) > and a{t) < 0, we have ('{t) ^ 1/e. We also see that C'(t) ^ whenever 
C{t) ^ l/2a_(_ and ('(t) ^ whenever ({t) ^ l/2a_. Since C(0) belongs to the interval 
[1/20+, l/2a_], C{t) must remain in this interval for all t. □ 

As we have already seen in ( p. 14 ), the probability of finding outside a strip of width 
much larger than ^2v{t) is very small. By Lemma 3.1, we now know that y^2u(t) behaves 
approximately like a\a{t)\~^/'^ . One of the key points of the present work is to show that 
the whole path {ys}o<s<t remains in a strip of similar width with high probability. The 
strip will be defined with the help of C{t) instead of v{t), because we need the width to be 
bounded away from zero, even for small t. 

To investigate we need to estimate the stochastic integral from ( 3.12 ). Lemma kA 
in the appendix provides the estimate 



P| sup / ip{u) dWu ^6} exp 



2 /o ^{u)^ dn 



(3.21) 
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for Borel-measurable deterministic functions ip{u). Unfortunately, this estimate cannot be 
applied directly, because in ( 3.12| ), the integrand depends explicitly on the upper integra- 
tion limit. This is why we introduce a partition of the interval [0,t]. 

Lemma 3.2. Let p : / ^ M+ be a measurable, strictly positive function. Fix K €N , and 
let = uq ^ ui < ■ ■ ■ < UK = t be a partition of the interval [0, t] . Then 



(3.22) 



k=l 



where 



Pk = e^p< --^{ inf p(s) e 



2 \«fc_l<S^Mfc 



2 2aiuk,s)/e\ i 



1 ruk 



2a(uk,s)/e 



-1 



Proof: We have 



tO,0 



sup -r^ 



> h 



} 



(3.23) 



(3.24) 



sup r e°(^'")/" dTy„ 

^o^s^tp{s) Jo 



a J 



^2^p0.o| sup / e-=(")/^dW„^ 



^a{s,u)/e 



inf p(s)e-"(^)/4. 



Applying Lemma A.l| to the last expression, we obtain ( 3.22D . □ 

We are now ready to derive an upper bound for the probability that leaves a strip 
of appropriate width hp{s) before time t. Taking p{s) = \/C{s) will be a good choice since 
it leads to approximately constant Pk in (|3.22D . 



Proposition 3.3. There exists an r = r(o+,a_) such that 

„0| . .1/^2 



sup ^ ^h]^ C(t,.)exp{---(1 - re)}, 

•-O^s^t a/C('S' ^ 1-2(7^ J 



where 



Proof: Let 



C{t,e)J^ + 2. 



K 



2^2 



For /c = l,...,iir — 1, we define the partition times by the relation 



(3.25) 



(3.26) 



(3.27) 



(3.28) 
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which is possible since a{t) is continuous and decreasing. This definition impUes in partic- 
ular that a{uk,Uk-i) = —le^ and, therefore, Uk — Uk-i ^ 2e^/a_. Bounding the integral 



in (|3.23| ) by Ci'^k), we obtain 



P,^exp{-^^ inf ^e'-iu.,sys\ (3.29) 



We have e2=("'='")/^ ^ e"'''^ and 



C(5) - Ciuk) = - I C'iu) du ^ (3.30) 

J s ^ 

Since C(^A:) ^ l/2a+, this implies 

P,<e.p{-ii!(l-42±.)e-}, (3.31) 



and the result follows from Lemma 3.2. □ 



Remark 3.4. If we only assume that a is Borel-measurable with a{t) ^ — a_ for all t ^ I, 
we still have 

pO'Oj sup ^ V\/2^| ^ C(t,e)exp|-J:^e-44. (3.32) 



To prove this, we choose the same partition as before and bound the integral in ( ^.23|) by 
e/2a_. 



We now return to the nonlinear equation (|3.5| ), the solutions of which we want to 
compare to those of its linearization ( 3.11| ). To this end, we introduce the events 



Cltih) = {l^: \ys{uj)\ < h^/Ci^ys G [0,t]} (3.33) 
O0(/i) = [uj: |yO(a;)| < /i^/cRVs G [0,t]}. (3.34) 

Proposition |3.3| gives us an upper bound on the probability of the complement of Q^{h). 
The key point to control the nonlinear case is a relation between the sets Clt and (for 



slightly different values of h). This is done in Proposition 3.6 below. 



Notation 3.5. For two events $7i and CI2, we write Cli C if ^-almost all w G ili 
belong to ^2- 



Proposition 3.6. Let 7 = 2^2a+ Mja^L and assume that h < dy/a-/2 A 7 ^ . Th 



en 



Utih) c' + (3.35) 

0°(/i) c' !^t( [1 + 7/1] /i)- (3.36) 
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Proof: 

1. The difference Zg = ys — Vg satisfies 

dz 1 

^ = - [a{s)zs + + z„ s)] (3.37) 

with zq = P-a.s. Now, 

zs = - r e^(^'")/^ + z^, u) du, (3.38) 
£ Jo 

which imphes 

\zs\ < - r e"(^'")/^ !%„,«) I (3.39) 
£ JO 

for all s e [0,t]. 

2. Let us assume that lo £ ^t{h)- Then we have for all s £ [0,t] 

\y^iu;)\^h^)^^^^, (3.40) 
^ la- I 

and thus by ( ^^391) , 

|..(^)Ki re-(^'«)/^^d^. (3.41) 

The integral on the right-hand side can be estimated by ( ^.16 ), yielding 

- e=(^'")/" du ^ 2C2e(s) ^ — . (3.42) 
e Jo «- 

Therefore, 

M/i2 My/a^h 



\zs{u;)\ ^ ^ ^ ^\/CM, (3.43) 

which proves ( 3.35| ) because |ys(<^)| ^ |ys('^)| + ks('^)|- 
3. Let us now assume that u> S ^ti^)- Then we have lysCi^)! ^ d/2 for all s G [0, t] as 
in ( p. 401 ). For 5 = 'jh, we have (5 < 1 by assumption. We consider the first exit time 

r = inf{se [0,t]: \zs\ ^ 6hy/C(7)} £ [0,t]U{oo} (3.44) 

and the event 

A = n^n{uj: t{lo) < oo} . (3.45) 



li u) £ A, then for all s £ [0,r(a;)], we have |ys(ci;)| ^ (1 + 5)hy/ C{s) ^ d, and thus by 
d^) and (|1|), 

|.,(.)| < 1 r e=<-.../« MLti)!^! a„ < MLti)!^! < v^R. (3.40) 
e ./n 2a_ 2o 



However, by the definition of r, we have (a;)! = 5h^y ({t{uj)), which contradicts 
(Q) for s = t{uj). Therefore P{^} = 0, which implies that for almost all uj £ fi^, 
we have |-Zs('^)| < Shy^({s) for all s £ [0,t], and hence 



|y,(w)| < (1 + 6)h^/Cis) Vs £ [0, t] (3.47) 
for these cu, which proves ( |3.36D . □ 
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We close this subsection with a corohary which is Theorem restated in terms of 
the process yt- 

Corollary 3.7. There exist ho and Eq, depending only on f , such that for e < Eq and 
h < ho, 

pO'Oj sup ^= >h}^ at, e) exp\-\^ [l - 0(e) - Oih)]}. (3.48) 



Proof: By Proposition |3.6| and Proposition 

pO,or g,jp > ^ pO,0 1 > ^\ 

^C{t,e) exp{-^^(l-re)}, 

where h = (1 + 'yhi)hi, which implies 



(3.49) 



hi = —[^/l + i^ - 1] ^ h[l - ^h] (3.50) 
where we have used the relation vT+~2x ^ 1 + x — ^x"^. □ 
3.2 Unstable case 



We now consider a similar situation as in Section |3.1| , but with an unstable equilibrium, 
that is, we assume that a{t) ^ ao > for all t € I. Theorem 2.1 shows the existence of a 
particular solution x^^^ of the deterministic equation (|2.6|) such that la;^^''* — x*{t)\ ^ cie 
for all t S I. We are interested in the stochastic process yt = xt — xf'^'^ , which describes 
the deviation due to noise from this deterministic solution x"^^^ . It obeys the SDE 

dyt = - [a{t)yt + b{yt, t)] dt + ^ dWt, (3.51) 
e ye 

where 

m = aeit) = d,f{x^^\t) 
b{y, t) = Uy, t) = /(5?f * + y,t)- /(xf t) - a(t)y 



are the analogs of a and b defined in ( |3.6D . Taking e sufficiently small, we may assume 
that there exist constants oq, ai, d > 0, such that the following estimates hold for all t G / 
and all y such that \y\ ^ d: 

a{t) ^ -ao, \a'{t)\ ^ ai, \b{y,t)\ ^ My\ (3.53) 



The bound on |a'(t)| is a consequence of the analog of ( |3.10| ) together with the fact that 
\x^^' -x*{0)\ = 0{e). 

We first consider the linear equation 

dy° = -a{t)y^ dt + ^ dWf (3.54) 
e ye 
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Given the initial value yQ, the solution at time t is a Gaussian random variable with 
mean e"^*^/'' and variance 



vit) 



2 rt 



e Jo 



(3.55) 



where a(t, s) = a(«) du ^ ao(* — s) for t ^ s. The variance can be estimated with the 
help of the following lemma. 



Lemma 3.8. For < e < 2a^/ai, one has 



e Jo 



2a(0) 2a{t) 



[l + 0{e 



(3.56) 



Proof: By integration by parts, we obtain that 







2a(0) 



e f a'{s) 



2a{t) 2 Jq a(s)2 



which implies that 

ft 



1 _ £^ 

2al 



1 + -^ 

2ag 



2a(0) 2a(i) 

By our hypothesis on e, the first term in brackets is positive. 



^2a{t,s)/e (3.58) 



□ 



Unlike in the stable case, the variance grows exponentially fast (at least with e^"o*/^). 
If p ^ Ii/qI, we have 



P°'^8{ sup |y°| < p} ^ P°'^°{ly?l < P} 

.p-yOe^W/^ ^-x^/2v(t) 



dx ^ 



2p 



(3.59) 



which goes to zero as pa~^ e~°'^^^^'^ for t — > 00. In this estimate, however, we neglect all 
trajectories that leave the interval (— />, p) before t and come back. We will derive a more 
precise estimate for the general, nonlinear case by introducing a partition of [0,t]. 

The following proposition, which restates Theorem 2.5 in terms of yt, is the main result 
of this subsection. 

Proposition 3.9. There exist constants £q, Hq > such that for all h ^ a A Hq, all e ^ eo 
and for any given yo with |yo| \/2a(0) < h, we have 



tO,yo 



sup\ys\^/Ms) <h} !^^/^exp\-K^^}, (3.60) 



where k = § (l - 0{h) - 0{e)) . 



22 



Proof: 

1. Let K £ N and let = uq < ui < • • • < uk = t he any partition of the interval [0, t]. 
We define the events 

Ak = \uj: sup \ys\y^2a{s) < h\ 



Let gfc be a deterministic upper bound on = P"*'^"fe{Afc}, valid on Bk- Then we 
have by the Markov property 



pO'J^oj sup \ys\^/Ms)<h\ 

A:=0 

A'"2 K-l 



2. To define the partition, we set 

for some 7 G (0, f] to be chosen later, and 



(3.63) 



/j2 

a{uk+i-,Uk) = '^e^-, k = 0,...,K-2. (3.64) 

Since a{uk+i,Uk) ^ ao('Wfc+i - ^ifc)> we have n^+i - Uk ^ s-'^d ^^ing Taylor's 

formula, we find for all s G [uk, u^+i] and all A; = 0, . . . , X — 1 

/i^ ai a(s) /i^ ai 

1- — 327e ^ ^ 1 + — 327e, (3.65) 

cj"= Qq a[uk) Cg 

where a\ is the upper bound on |a'|, see (|3.53D . In order to estimate P^, we introduce 
linear approximations {y\ )te[Mfe,Ufc+i] for A; G {0, . . . , -fC — 2}, defined by 

dyf ) = ia(t)yf ) + ^ mf\ = (3.66) 

fc ye 

where Wl > = Wt - Wu^ is a Brownian motion with WuJ = which is independent 
of {Ws : ^ s ^ Uk}. If a; G A^^ we have for all s G [ufc, lifc+i] 



|y.(a;)-yW(u;)| ^ - f\''^'^-^l'\b{yu^u)\du 

^ — [l + 0{s)] ^ro- 



(3.67) 



2ao a('Ufc) y^2a(s)' 
where tq = Me(2aQ)~^/^ + 0(6). This shows that on Ak, 

|yf )(u;)| ^ [1 + roh] — i= Vs G K, u^+i]. (3.68) 
Y^2a(s) 
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3. We are now ready to estimate Pk- (|3.68|) shows that on B/., 



Pk ^ P"fe'f-fc I sup \yf^ 



W2a{s) < /i(l + ro/i)} 



p«.,j/.,||y(fc)^j^2a(n,+i) < Kl + roh)} 
1 2/i(l + ro/i) 



(3.69) 



where v^,}^-^ denotes the conditional variance of yn^'^^.i, given yu^. As in ( 3.56 ) 



e L,. 2 



a(wfc+i) 



[1 + Oie)]. (3.70) 



It follows that 



1 



a{uk) 

i + 27^)(i-ij;;57. 



[1 - 0(e)] 



[1 - 0(e)] 



(3.71) 



l-|2(l + 27)ej [l-O(^) 



^7^'[l-0(e)]- 

Inserting this into ( |3.69 ), we obtain for each k = 0, K — 2 on Bi^ the estimate 

1 



2tt 727/? 



[l + 0{e)+0{h)]=:q. (3.72) 



Note that for any 7 G (I/tt, 1], there exist ho > and eo > such that q < I for all 
h ^ ho and all e ^ eq. Since qk-i = 1 is an obvious bound, we obtain from (3.62) 



sup \ys\^/m7)<h} ^ g^-i ^ -expj- 



g2log(l/g2)|. (3.73) 



Choosing 7 so that q'^ = 1/ e holds, yields almost the optimal exponent, and we obtain 



pO'^oj sup \ys\VM^ <h}^ V^exp{-K^^}. 



(3.74) 
□ 



4 Pitchfork bifurcation 

4.1 Preliminaries 

We consider the nonlinear SDE 



dxt = -f{xt,t)dt + ^dWt 



in the region M. = {(a;,t) € M^: |x| ^ d, \t\ ^ T}. We assume that 



(4.1) 
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• there exists a constant M > such that /(x, t) is three times continuously differen- 
tiable with respect to x and t and \dxxxf{x,t)\ ^ 6M for all {x,t) G M.; 

• = -f{-x,t) for all G 

• / exhibits a supercritical pitchfork bifurcation at the origin, that is (after rescaling), 

5a./(0,0)=0, 9t,/(0,0) = l and a,,,/(0, 0) = -6 (4.2) 

Using Taylor series and the symmetry assumptions, we may write for all (x, t) £ Ai 

f{x,t) = a{t)x + b{x,t) = x[a{t) + go{x,t)] 
dxf{x,t) = a{t) +gi{x,t) 

where a(t), go{x,t), gi{x,t) are twice continuously differentiable functions satisfying 

a{t) = dxf{0,t) = t + Oit^) 
go{x,t)=[-l + -fo{x,t)]x^ \go{x,t)\ Mx^ (4.4) 

gi{x,t)= [-3 + 7i(x,t)]x2 \gi{x,t)\ <3Mx2, 

with 7o,7i some continuous functions such that 70 (0,0) = 71(0, 0) = 0. The following 
standard result from bifurcation theory is easily obtained by applying the implicit function 
theorem, see |GH, p. 150] or [|^, Section II. 4] for instance. We state it without proof. 

Proposition 4.1. If T and d are sufficiently small, there exist twice continuously differ- 
entiable functions x*, x : (0, T] — > R + of the form 

X*{t) = Vt[l + OT{l)\ 

r— (4.5) 

x(t) = 7t73[l + OT(l)] 

with the following properties: 

• the only solutions of f{x,t) = in Ai are either of the form (0,t), or of the form 
{±x*{t),t) with t > 0; 

• the only solutions of dxf{x,t) = in M are of the form (ibx(t),i) with t ^ 0; 

• the derivative of f at ±x*(t) is 

a\t) = dxf{x\t),t) = -2t[l + Ot(1)] . (4.6) 

• the derivatives ofx*{t) and x{t) satisfy 

^ = -L[l + Ot{1)], ^ = -^[1+Ot{1)]. (4.7) 
at 2Vt at 2V3t 



As already pointed out in Section 2^, there is no restriction in assuming T and d to 
be small. Thus we may assume that the terms C't(I) are sufficiently small to do no harm. 
For instance, we may and will always assume that a*{t) < 0. 

Equation ( |4.4| ) also implies the existence of constants a+ ^ a_ > such that 

a+t ^ a(t) ^ a^t for -T ^ t ^ 

(4.8) 

a^t ^ a{t) ^ a+t for < t ^ T. 
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The function a{t, s) = a{u) du thus satisfies 



-ia+(s^ -t^) ^ ait,s) ^ -ia_(s^ -t^) 
ja_t^ — ia+s^ ^ a{t,s) ^ ia+t^ 



if s ^ t ^ 
if s ^ ^ t 
if ^ s < t. 



(4.9) 



We are going to analyse the dynamics in three different regions of the (t, x)-plane: near 
X = for t ^ y^, near x = for t ^ y^, and near x = x*{t) for t ^ -ye. In order to 
dehmit the last two regions, we introduce (somewhat arbitrarily) the function 



set 



and define the region 



x{t) = \/Ax*(i), 



a{t) = dj{x{t),t), 



V = {(x,t): y/£i^ti^T,\x\< x(t)}, 



which has the following properties: 

(a) for all (x, t) ^ D with x ^ 0, one has 

-f{x,t) ^ Ka{t) 



(4.10) 



(4.11) 



(4.12) 



with K = 1 — A — C't(I) 



(4.13) 



(b) for all (x, t) G [-d, d] x T] \ V, 



dxf{x, t) ^ a{t) ^ -r/a(t) with = 3A - 1 - cir(l)- 



(4.14) 



For our results to be of interest, k > and ?7 > are necessary, which requires A € (|, 1). 
As we shall see, we will actually need A € (^, \). Furthermore, in Section we need to 
assume that cr | log cr|^/^ = 0{y/e). 



In the following subsections, we investigate the three different regimes: In Section 4.2 , 
we analyse the behaviour for t ^ ^/£. Theorem is proved in the same way as Theo- 
rem |2.3| , the main difference lying in the behaviour of the variance which is investigated 
in Lemma 14. 2|. 



Section ^4.3| is devoted to the rather involved proof of Theorem p.Qj We start by giving 
some preparatory results. Proposition 4.7 estimates the probability of remaining in a 



smaller strip 5 in a similar way as Proposition 3.£. We then show in Lemma 4.8 that the 



paths are likely to leave T> as well, unless the solution of a suitably chosen linear SDE 
returns to zero. The probability of such a return to zero is studied in Lemma [4.9|. Finally, 



Theorem 2^ is proved, the proof being based on an iterative scheme. 

The last subsection analyses the motion after tt>- Here, the main difficulty is to 
control the behaviour of the deterministic solutions, which are shown to approach x*(t). 

We then prove that the paths of the random process are likely 

The proof is similar to the 



cf. Proposition 4.11. 
to stay in a neighbourhood of the deterministic solutions. 



corresponding proof in Section 3.1 
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4.2 The behaviour for t ^ ^/e 

We first consider the linear equation 

dx^t = -a{t)x^ dt + ^ dWt 
£ ye 

with initial condition x^^ = xq at time to £ [~^)0). Let 



,{t,to) = — re2a(t,s)/e^^_ 

e Jto 



(4.15) 



(4.16) 



denote the variance of x^. As before, we now introduce a function C(t) which will allow us 
to define a strip that the process xt is unlikely to leave before time ^/£, see Corollary 4.5 
below. Let 



,2a{t,to)/e 



+ - 



1 



2\a{to)\ £ 
The following lemma describes the behaviour of C,{t) 



(4.17) 



to 



Lemma 4.2. Assuming e ^ 4a(to)^ A (to/2)^, there exist constants c± = c-i-(a+,a_) such 
that 



for to^t ^ - 
/or -Ve ^t^^/e 
for^/£^^t^^T. 



(4.18) 



moreover, a'{t) > on [to,t], then C(t) is increasing on [to,t]. 

Proof: The upper bounds are easy to obtain. For to ^ ^ ^ have, using t^-s^ ^ 

2t{t - s), 



to 



2|a(to)| \t\ L2a_ 2a^ 



For — \/e ^ t ^ 0, the hypothesis e ^ 4a(to) implies 



^(t) ^ _e-"- / e-'^-^'/^ds + 



1 



1 



to 



2\a{to)\ 



dn + 1 



(4.19) 



(4.20) 



For ^ t ^ ^/e, a similar estimate is obtained by splitting the integrals for s ^ and 
s ^ 0. For t ^ ^/e, we have 



-2a{t)/e 



du + 



du + 1 



(4.21) 



To obtain the lower bound, we first consider the interval to ^ t ^ ^to, where we use the 
estimate t^ — ^ 2to(t — s), valid for all s £ [to,t], which yields 



^-2a+\to\(t~s)/e ^ ^ _ ^ 



1 



to 



2a+\to\ 



2a+\t\ 



(4.22) 
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For ^ i ^ — we have — 3t{t — s) for all s E [2t, t], and thus 



1 - e-3"+ 
3a+|i| 



(4.23) 



where we used the relation Iq ^ in the last step. By the same relation, we obtain 



-2a{t)/e ^^^^ ^ 



1 

1 



-1 



e-''+" du 



e~"+" dn 



for - i^t^:^/e, (4.24) 
for t ^ ^/e. (4.25) 



Finally, assume that a'{t) > for all t, and recall that C(0 is the solution of the initial 
value problem 



dC_2a(t) 1 1 

d^-^^+e' ^^*°^-2K^ 



(4.26) 



Since C(t) ^ 0, C' > for all positive t. For negative t, C,' is positive whenever the function 
V{t) = C{t) + l/2a(t) is negative. We have V{to) = and 

dV_2a(tl 

dt ~ e ^ 2a(t)2- ^^■^'> 

Since V < whenever V = 0, V can never become positive. This implies C' ^ 0. □ 

The following proposition shows that the solution of the linearized equation ( [4.15 ) 
is likely to track the solution of the corresponding deterministic equation. 

Proposition 4.3. Assume that —T ^ to < t ^ yje. For sufficiently small e, 

sup ^^^7^= ^>a^C(t,e)exp{-^^[l-r(e)]}, (4.28) 



where 



C(t e) = 1"^^'^°^! + «+ + + 4 



and w/iere r(e) = ©(e) /or to ^ ^ ^ —\/^, O'^d r(e) = ©(a/e) for — a/e ^ t ^ a/e- 



(4.29) 



Proof: Let to = ""o < ■ ■ ■ < uk = t be a partition of the interval [tojt]- By Lemma 3.2 , 
the probability in ( [4.28 ) is bounded by 2^^^ P^, where 



P. = exp|---^ — 



C{Uk) Uk 

If t ^ —\fE^ we define the partition by 
-a(t,to) 



inf C(^^)e^"("^'")/^|. 



K 



2e2 



a{uk,to) = 2e k for k = 0, . . . , K — 1. 



(4.30) 



(4.31) 
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Estimating as in the proof of Proposition we obtain 

1 2e 



r i /I" / 
P.^exp|---^(1 



4£ 



}■ 



(4.32) 



Therefore, (|^) holds with C{t,e) = \a(t,to)\/e'^ + 2. 

For ^ t ^ y^, we define the partition separately in two different regions. Let 



-a{-^/e, to) 
2e2 



K = Ko + 



(4.33) 



The partition times are defined via 
-a{uk,to) = 2£^k 



for ^ A; ^ i^o - 1 
iov K^i^k i^K -I. 



(4.34) 



In the first case, we immediately obtain the bound ( 4.32| ). In the second case, estimating 
Pfc in the usual way shows that 



Pk ^ exp{-i J(l - ^[l + 2a+c+]) e-^+^j. 

Finally, let us note that, for — -y/e ^ t ^ -v/e, 

,^^ K-x^,to)| ^2 k^^^^^^ 

which concludes the proof of the proposition. 

Let us now compare solutions of the two SDEs 



(4.35) 

(4.36) 
□ 

(4.37) 
(4.38) 

(4.39) 
(4.40) 



Proposition |4.3| gives us an upper bound on the probability of the complement of il^(/i). 
We now give relations between these events. 

Proposition 4.4. Let t G [to,\/e] ^'^^ l^ol ^ hje^l'^, where we assume h? < e/7 for 



Ax\ = -a(t)x? dt + ^ dWt 
£ ye 

dxt = -f{xt,t)dt + ^dWt 
£ ye 



Xto = Xo, 



where to ^ We define the events 

n^h) = {w: |2;°(u;) -xoe"("'*«)/"| ^ h^/C{sjys G [to,t]^ 
nt{h) = [uj: |x,(a;)-Xoe°(^'*°)/^| <^hy^yse [to,*]}- 



7 = M(l + 2^fc+/^ and ^ d^^l{\ + . Then 



nt{h) c n'l 



1 + 7- 



1 + 7- 



e 



h 
h]. 



(4.41) 
(4.42) 
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Proof: Assume first that lo G ^^'^ l^t 6 = ■yh'^/e. Then we have 5 < 1 by 

assumption. By (^^), the difference „ satisfies 



to 



We consider the first exit time 



r = inf{sG [to,t]: \zs\ 6h^/C{s)} £ [to, t] U {cx)}. 

For all u! in the set 

A = n^{h) n {uj: t{uj) < oo}, 
and s G [tQ,T{uj)], we have by the hypotheses on h and xq together with Lemma 



,H| ^ |xo|+/iVC(s) ^ (l + (l + ^)^);T7i ^rf- 



Therefore, yields 



+ + 



/4 



3 1 



(4.43) 



(4.44) 



(4.45) 



(4.46) 



to 



The integral is bounded by 2(2e{s), which can be estimated by Lemma 4.2 
Thereby, we obtain 



c+ 



(4.47) 
once again. 

(4.48) 



^ e 

which leads to a contradiction for s = t{uj). We conclude that P(A) = 0, and thus 
t{uj) = oo for P-almost all u) G ^l^{h). This shows that |2;s(w)| < 6hy^^(s) and thus 
\xs{uj) — e"^*'*'')/^! < (1 + 5)h\/ C(s) for all these uj and all s G [to;^]) which proves 
( |4.42[ ). The proof of the inclusion ( ^.41| ) is straightforward, using the same estimates. □ 

The two preceding propositions immediately imply the main result on the behaviour 
of the solution of the nonlinear equation ( [4.38 ) for t ^ y/e, i.e.. Theorem 2^, which we 
restate here with an arbitrary initial time to £ 

Corollary 4.5. Assume that —T ^ to < t ^ ^fe. Then there exists an Hq > such that 
for all h ^ ho^/e and all initial conditions xq with \xq\ ^ hje^l^ , the following estimate 
holds: 



sup 

to^s^t 



\Xs - xoe' 



a(s,to)/e| 



>h]^ C(t,e)exp{-i^[l - r(e) - 0{hye)]], 



(4.49) 



where C{t,£) and r(e) are given in Proposition U 
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4.3 Escape from the origin 

We now consider the SDE (14.11), written in the form 



Axt = - [a(t)xt + h{xt,t)\ dt + ^ dWi, (4.50) 
£ ye 

for t ^ to ^ y^, where we assume that {xtgl ^ x{to). Our aim is to estimate the first exit 
time Tx> of xt from T? defined in ( |4.12[) . We recall that a{t) + ^b{x,t) ^ Ka{t) in T), see 
( 14.131 ). Moreover, we have a_t ^ a{t) ^ a+t, ^ a'{t) ^ ai, and |6(x,t)| ^ M\x\^ i n P. 
We first state a result allowing to estimate the variance of the linearization of (l4.50| ) . 

Lemma 4.6. Let a{t) he any continuously differentiahle, strictly positive, increasing func- 
tion, and set a{t,s) = f^a{u)du. Then the integral 

fj2 ft 

v{t,s) = — e2"(*'")/^dn (4.51) 

^ J s 

satisfies the inequalities 

^[eM...,/._i]^„(,.,^_^3..,.,.,A. (,.5,) 

Proof: Using integration by parts, we have 

^ ^ l2a{s) 2a{t) J, 2a{uf J ^ ' 

The upper bound follows immediately, and the lower bound is obtained by bounding the 
exponential in the last integral by 1. □ 

Our first step towards estimating rx) is to estimate the first exit time ts from a smaller 
strip 5, defined as 



S=[{x,t): V^^t^r,|a;| < 

I Va(s)J 



(4.54) 



where we will choose 



h = 2o-V|loga|. (4.55) 



Proposition 4.7. Let to ^ \/e dnd |a;o| ^ h/ y^a{to)- Then, for any /x > 0, we have 
under the condition 
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Proof: 

1. For K £ N , we introduce a partition to = uq < ■ ■ ■ < uk = t oi the interval 
[to,i], which will be chosen later, and for each k, we define a linear approximation 

y^t Jte[«fc,ufe+i] Dy 

dxf ) = \a{t)xf^ dt + ^ dTyf ) 4? = X.., (4.58) 



Lemma 4.6 



where = Wt — Wu^.- Assume that \xs\\/a(s) ^ h for all s G [u^jftfe+i]. Then by 



1 (4.59) 
for s £ [uk,Uf:^i], provided the partition is chosen in such a way that for all k 



MVa(nfe+i) ° 



2. If \xuj\/a{uk) ^ /i, then we have 

P"'='""4 sup ^ /i| ^ P"'^'""^ I I VaCufe+i) ^ 2h} 



4h 

v2^^V+i«(^fe+i) 



where the variance 

^2 l'Uk+ 



W =_ / e2aK+i,.)/e^^ ^4g2) 



can be estimated by Lemma 4.£. We thus have by the Markov property 

P = P*0'^n| sup |2;,|v^^4^ n ( ^ ^ ^1 V (4.63) 

3. We now choose the Uk in such a way that fi^^^^a(nA:+i) is approximately constant. 
Given > 0, let 

i=h^(^Y (4.64) 

(Observe that ^ ^ S/i^/vr > cr^/2.) Choosing iiT as the smallest integer satisfying 

2a(t,to) 

^ ^ slog(2£/a^) ' ^'-''^ 
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we define the partition by the relations 

a{uk+i,Uk) = ^log^, 
< a{uK,UK-i) ^ -log^. 



for A;e{0, ...,K-2}, 



Then we have 
4 h 



^CT^2£/cj2 -1 



exp 



«(Mo) log[(g)"-f5a 



(4.66) 
(4.67) 

(4.68) 



which proves ( [4.56| ). 



4. It remains to show that condition ( 4.60 ) is satisfied. Since 



a{uk+i) ai ^ ^ 1 , 1 

^ 1 + — — -{uk+i - t/fc) ^ 1 + 71-75^ log 



a{uk) a{uk) 
the condition reduces to 



+ 



log 



16 //l^Xl+M 



vr 



16 



M 4 0-2^ 



(4.69) 

(4.70) 
□ 



which is satisfied whenever condition ( [4.57| ) is satisfied. 

We want to choose n in such a way that P*o,a;o{^^ ^ t} < (/i/a)^ e"'""^*'*")/^ holds 
with the same k as in ( 4.13 ). We opt for fi = 2, because this choice guarantees the 
above estimate for all possible k without choosing a K-dependent fi. For h = 2(Ty^|log a\, 
Condition ( 4.57] ) becomes a consequence of the following slightly stronger condition 

a|log(T|3/2 



(4.71) 



which we will assume to be satisfied from now on for the rest of this subsection. 

The second step is to control the probability that xt returns to zero after it has left 
the strip S. To do so, we will compare solutions of ( 4.50 ) with those of the linear equation 

(4.72) 



dxj = -ao{t)x^ dt + dWt, 
£ ye 



where ao(t) = na{t) satisfies ao(t) ^ f{x,t)/x in D. The following lemma shows that this 
choice of ao(s) implies that \xs\ ^ \x^\ holds as long as Xg does not return to zero (Fig. ^). 
This implies that if x^ does not return to zero before time i, then Xg is likely to leave V 
before time t without returning to zero. 

Lemma 4.8. Let to ^ O'nd assume that < xq < x{tQ). We define 

p+(t) = {{x,s): ^/ei^sii,tandQ<x< i{s)] (4.73) 

and denote by r-p+ the first exit time of Xs from V^it). Let he the time of first return 
to zero of x^ in [to,t], where we set = oo if x'^ > for all t G [to,t]. Then Xg ^ for 
all s ^ r-p+ A t and 

P*O'^o|0 < X, < x(s) VsG [to,t],T° = oo} ^ P*«'^«|0 < x° < x(s) Vs E [to,t]} 

^-Ka(tM>)/e (4.74) 



x{t)\/aQ{t) 



I -no 
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Figure 3. Assume the path xt exits the region S at time ts, say by passing through the 
upper boundary of S. We introd uce a process Xj, starting on the same boundary at time 
Tg, which obeys the hnear SDE ( 4.72| ). Let r° be the time of first return to zero of x^. 



Then xt lies above x^ for Tg < t ^ t° . In case Xt also becomes negative, the two processes 
may cross each other. The probability of x^ ever returning to zero is bounded by ct^*^. If 
x° does not return to zero, Xt is likely to leave V. 

Proof: 

1. Let g{x, s) = f{x, s) — ao{s)x. By assumption, g{x, s) is non-negative for (x, s) S V'^. 
The difference g satisfies the equation 

1 

Zs = zto + - / [g{xu,u) + ao{u)zu\ du (4.75) 

with ztf, = 0. Since g{xs, s) ^ for to ^ ^ r-p+ A t, 

1 P 

^ ^;to + - / ao{u)zudu, (4.76) 



Jo 

follows for all such s and, therefore, Gronwall's lemma yields 

Zs ^ zto e'^°("'*»)/^ = for all s G [to, Tjy+ A t]. (4.77) 

This shows Xg ^ for those s. Now assume Tp+ = 00 and r'' = 00. Then, (|]73) 
implies that < ^ < x(s) for all s ^ t, which shows the first inequality in 



(4.74). 



2. Xg being distributed according to a normal law, we have 



{0 < x° < x(s) Vs G [to,t]} ^ P*0'^'o{0 < x° < x(t)} 



^ x(t) (4.78) 
where the variance vo(t,to) can be estimated by Lemma [4.6|. This proves the second 



inequality in ( |4.74D . □ 



The previous lemma is useful only if we can control the probability that the solution Xj 
of the linearized equation returns to zero. The following result estimates this probability 
and its density. 



34 



Lemma 4.9. Let to ^ y/e and assume that x^^ = P > o"/y^ao(to)- Denote by the time 
of the first return of to zero. Then we have 

-P*0'^{rO<t}^— 7^^e-»«(*°)^/'^ -Vao(t)ao(to)^ ' ^ ^^^^ ■ (4.80) 
Proof: 

1. Since by symmetry, P'^°'''{xj ^ 0} = ^ on {r'' < t}, we have by the strong Markov 
property 

F'°'P{x^t ^ 0|r° <t} = l- (4.81) 

We now observe that 

F^°^P{x°t ^ 0} = F^'>^P{x°t ^ 0, r° ^ t} + P*«'^{x? ^ 0, t° < t} 

= P*o,P{t-0 ^ t} +p*o,P{^0 ^ Q|^0 ^ ^|pto,P|^0 ^ ^1 

= 1 _ pto,P{^0 ^ ^1 ^ lpto,P{^0 ^ ^1 (4-^2) 

= 1 - iP*0'''{T° < t}, 

which imphes 

P*0'P{r° <t} = 2[l- F^°'P{x^ ^ 0}] = 2F*°'P{x^ < 0}. (4.83) 

2. Next, we use that x^ is a Gaussian random variable with mean pe'^"^*'*"^/^ and variance 

^2 rt 

y^(t,to) = — e2«"(*'*)/^ ds. (4.84) 

£ Jtn 



By Lemma 4.6 



2 2K.ait,to)/e 2 



and we thus have 



1 r° r - «e'*"(*'*o)/^'l^ -1 
P*o,P{^o < o| ^ ^ / expl-^^^^-^^ ^} dx 

^2TTVoit,to) J-oc ^ 2vo{t,to) J 



^j,Ka(t,to)/e 



27r J-oo 2 



e-^/My^-e-=, (4.86) 



which proves ( 4.79| ), using (|4.83| ) and (|4.85|) 



3. In order to compute the derivative of F^^'P{x^ < 0}, we first note that 

^^vo{t, to) = Y + ^^^Mt, to). (4.87) 
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Differentiating tlie second line of ( [4.86 ), we get 



dt 



P*0'''{j;° < 0} 



1 



27r 
1 



exp 



p2 ^2K.a{t,to)/e 
^Ka{t,to)/e 



dt 



pe' 



Ka(t,to)/e 



\Ao(Mo) 



27r 2 e vo{t,tof/'^ 

1 1 cr^ Q-Ka{t,to)/e _ 



■r, e 



2ttP e y/voit,to) 
27r cr e 



(4i 



-2Ka{t,to)/e 



-2Ka{t,to)/£ 



where we have used the facts that H > ao(to)p^/<7^ > 1 and that He ^ is decreasing 
for H > 1. Now, foUows from (|^). □ 



Assume for the moment that starts "on the border" of S, i.e. in p(to) = h/ ^Ja{tQ) = 
y/nh / yaoito] . Then, by our choice h = 2ay|loga|, Estimate ( [4.79| ) shows that the 
probabihty for to return to zero cannot exceed e~'^'^^*°^^ ^'^ = a^'^. 

We are now ready to prove the main estimate on the first exit time rp, which is the 
most important of our results. Since the proof is rather involved, we restate Theorem 2.E 
here for convenience. 



Proposition 4.10 (Theorem 2.S). Let to ^ y/e and \xq\ ^ x{tQ). Then 



|logcj| 



1 + 



a(t,to] 



-Ka{t,tci)/e 



y^l — Q-2Ka{t,to)/£ ' 



(4.89) 



where Cq > is a (numerical) constant. 



The strategy of the proof can be summarized as follows. The paths are likely to leave 
5 after a short time. Then there are two possibilities. Either the solution of the linear 



equation ( [4.72|) does not return to zero, and Lemma 4.8 shows that xt is likely to leave D 
as well. Or x^ does return to zero. Using the (strong) Markov property and integrating 
over the distribution of the time of such a (first) return to zero, we obtain an integral 
equation for an upper bound on the probability of remaining in D. Finally, this integral 
equation is solved by iterations. 



Proof of Proposition 4.10 



1. We first introduce some notations. Let 

<^>tis,x)=F'^^Tv^t} 



sup — 



<1 



(4.90) 



and define p{t) = h/ \J a{t). We may assume that p{t) ^ x{t) for all t (otherwise we 
replace x by its maximum with p). For t ^ s ^ ^fe we define the quantities 



qt{s)= sup $t(s,x), 

\x\^p{s) 

Qtis) = sup $t(s,x). 

p{s)^\x\^x{s) 



(4.91) 
(4.92) 
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2. Let us first consider the case ^ p{s)- Recall that S = {{x,t): \x\ < p{t)}. By 
Proposition and the strong Markov property, we have the estimate 

=P^'"{r5 ^t}+P^'"|r5 <t, sup ^<l| 

I rs<^u<itXiu] ) 



^ (^!)%— (M)/.+]E^,-|lj^^^^(^^)Q^(^^)}. (4.93) 



The second term can be estimated by integration by parts, see Lemma A. 2, Let 
Qtiu) be any upper bound on Qt{u) satisfying the hypotheses on g in that lemma. 
Since Qt{u) ^ Qt{t) = 1, we may assume that Qt{t) = 1. Application of (A/7) with 
G{u) = 1 — (/i/o")^ e"'^"^"''')/'' shows that the second term in ( [4.93 ) is bounded by 

We have thus obtained the inequality 

q^^s) ^ 2(^)%-'^-(*'^)/^+At(^)'^*a(7x)^e-''"("'^)/^d«. (4.95) 

3. Consider now the case \x\ S [p(s),x(s)]. Since x i-^ fix,t) is an odd function, 
^t{s,x) = ^tis,—x) follows. Hence we may assume that x > 0. We consider the 
linear SDE (|4.72[ ) with initial condition x'^ = x, and denote by the time of the first 
return of x^ to zero. Then we have 

Ms,x) = F''^\t° ^ t, sup ^ < l| +P^'^|rO < t, sup ^ < l|, (4.96) 
1^ .<,<,,<t xiu) J I- xiu] ) 



and Lemma 4.8 yields 

r...|,o ^ M < iW ^")y^ , (4.97) 

>■ si^u^t X[U) ) ^T^a _ Q-2Ka{t,s)/e 



The second term in ( [4. 961 ) can be estimated using the density of the random variable 
r^, for which Lemma |4.9| gives the bound 

d n 1 2k3/2 /i ,2, 2 a(M) ^-'^'^'^M/e 

^^o{u) = — P^'^|r° < u| ^ — e"'^^ ^ , (4.98) 

du ^ -y/vr a e ^\ _ g-2Ko(M,s)/e 

We obtain 

P^'-{rO<i, sup ^^^|^jE.,.ri p.o,,^o| ^<lU 

s^u^t X{U) J I ^ ^ ^T^^u<iA X{U) J J 

V'ro(^t)^i(^i, 2;„)du 
^ / Vro(n)[gt(tx) +Qt(n)] du. (4.99) 
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4. Before inserting the estimate ( [4.95| ) for qt{u), we shall introduce some notations and 
provide bounds for certain integrals needed in the sequel. Let 



9(t,s) 



_ Q—2Ka{t,s)/e 

and (p = e-'="(*'^)/^ Then 

^ e— ('^'^)/^ g{u, s) dn ^ f ^9{u, s) dn 



(4.100) 



e 



e 



vr 2 

^ — ^ - 
2k k 



dx 



vr , 2 
< 



2k 7o v^x(1 - x) 2k k 



(4.101) 
(4.102) 



a[u] 



1-X2 



■ dx 



^ - log - ^ 



1 ^ a(t,g) 



01 i + v/T^ 

log —;^=^= 

e-'^"(i.^)^^ (4.103) 



where we used the changes of variables e 2Ka(«,s)/e _ ^(^^ — 0^) + ^ in ( 4.102| ) and 



X 



^ _ g-2«a(t,«)/e (|4l03D. 



5. Now we are ready to return to our estimate on ilj^o{u)qt{u) du, compare ( [4.99 ). 
Inserting the bound ( |4.95 ) on qt{u) yields two summands, the first one being 



V',o(n)e-'^"(*'")/" du 
4^3/2/ /in 3 _^/^2/^2 /■*a(n) e-2i^»M/s 



£ Y^]_ g— 2Ko(n,s)/e 



^2^f-f e-'^'^'/'^%-'^"(*'^)/^ 



(4.104) 



where we used ( [1.101| ) to bound the integral. The second summand is 

^%,o(^x) rQ,(^)^e-''"(^'")/^d7;d^z 

a(?;) 



K| - 



ye-^'^'/-' / g,(„)^e-"('''^)/^d^;, 



(4.105) 



where we used ( [1.101| ) again. 

We can now collect terms. Introducing the abbreviations 



rx(t)v^^ . and c 
the previous inequalities imply that 
Qt(s)^C5(t,s)+ce-''"(*'^)/^ 



/l\3 

TTKl —J e 



(4.106) 



[l + 5(u,s)] du. (4.107) 
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6. We will now iterate the bounds on Qtis). This will show the existence of two series 
{flnln^i and {bn}n^i such that 



Qtis) ^ Cg{t, s) + an e'''^^*'^)/^ +6„ Vn. 
To do so, we need to assume that 



e K 



; air, to) , 2 



By our choice ( [4.55| ) of h, this condition reduces to 

a2«|loga|3/4 = 0(VF), 



(4.108) 



(4.109) 



(4.110) 



which is satisfied for small enough e by our assumption ( 4.71| ) on cr, provided k > 1/2. 
Using the trivial bound Qt{u) = 1 in (f4.107| ), we find that ( |4.108| ) holds with ai = c 
and bi = cic/K. Inserting ( 4.108D into ([4.107D again, we get 



Qt{s) ^ Cg{t,s)+ce 



-Ka{t,s)/e 



t r- 



+ C 

< Cg{t,s) +c 
By induction, we find 



Cg{t,u)+an e-''"(*'")/= +6n 



a{u) 



-Ka{u,s)/e 



[l + g{u,s)] du 



l + C7(^ + ^UaJ^ + ^ 



3c. 



ttn+l = C 



1 + C 



1 + C 



a{t, s) 

£ 

a{t, s) 



n-l 



dm 



j=Q 



£ K 



+ C 



,a(M)^2 



+ - 

£ K 



(4.111) 
(4.112) 



as a possible choice, where we have used the fact that c(a(t, s)/e + 2/ n) ^ ^ by the 
hypothesis ( 4.10S| ). Taking the limit n ^ oo, and using c ^ -I ^ we obtain 

Qt{s) ^ Cgit, + + 3C) e-'^"(*'^)/^ ^ SCgit, s). (4.113) 

In order to obtain also a bound on qt{s), we insert the above bound on Qt{s) into 
( 4.95| ), which yields 



q^{s) ^ 2(^)%-'^"(*'^)/^ -^^"^{^y r^e-''"("'^)/^g(t,n)dn 



^ , 1 a(t, s) 
2 + 3kC - + ^^ 



/l\2 



-Ka{t,s)/e 



(4.114) 



by ( 4.103| ). This proves the proposition, and therefore Theorem 2.£, by taking the 
sum of the above estimates on qtis) and Qt{s). □ 
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4.4 Approach to x*{t) 

We finally turn to the behaviour after the time t = tx> > ^/e, when xt leaves the set P. 
By symmetry, we can restrict the analysis to the case Xr = x{t). Our aim is to prove that 
with high probability, xt soon reaches a neighbourhood of x*{t). 

We start by analysing the solution x^'^ of the deterministic equation 

with initial condition Xr''*'^ = x{t). 
Proposition 4.11. For sufficiently small e and T, 



e 



+ (x*(r) - x(r)) e-^"(*'")/" 



^ x*(t) - xP'" ^ C 
^ - xf ^ (x^^*'^ - x(r)) e-''"^*'^)/^ 



(4.115) 



(4.116) 
(4.117) 
(4.118) 



for all t G [t, T] and all r G [y^, T], where C > is a constant depending only on f. 



Proof: 



1. Whenever x^'^*'^ = x*{t), we have 



d / . dpi T\ dx*(t) „, ^, X X dx*(t) 



dt 



dt 



dt 



(4.119) 



which shows that x^°*'^ can never become larger than x*(t). Similarly, whenever 



det,' 



X, 



x{t), we get 



s^^{xr''^-m=mt)^t)-e^ 



(1 - A)i=^/2 ^ ^^(;l)] - [1 + Ot(1)] > 



(4.120) 



provided A < ^[1 — C't(I)]) which shows that x^^*'^ can never become smaller than 



':{t). This completes the proof of ( 4.116| ). 



2. We now introduce the difference = x*{t) — x'^^^''^ . Using Taylor's formula, one 



det,T 



immediately obtains that yf*^*'^ satisfies the ODE 

dy 



where 



e-^ = a*it)y + b''{y,t) + ex'''it) 



a*{t) ^ -alt 

^ b*{y,t) ^ M*Vty'^ 
K* 

."it) < ^, 



(4.121) 



(4.122) 
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with = 2[1 + Ot(1)], M* = 3[1 + Ot(1)] and K* = ^[1 + Ot(1)]. We first consider 

the particular solution yf*^^ of ( |4.12lD starting at time 4:^/£ in = 0. By (4.119), 
we know that yf'^^ ^ for all t ^ 4:^/e. We will use the fact that 

\fs Jt vS 

^ V . . 

4e /"^ e-" . e ^ ' 



where cq = S/og. We have used the transformation s = t — 4:£u/{aQt), introduced 
^ = aQt^/4e and bounded the last integral by 2. We now introduce the first exit time 
f = inf{t ^ A^/e: yf''^ ^ coet"^/^}. For 4^ ^ t ^ f , we have 

a\t)y + b\y,t)i^ [-a^t + Vtco^)y -al{l - "^^Yy. (4.124) 

Since M*/(ao)^ = |[1 + c(l)], the term in brackets can be assumed to be larger than 
i. Hence ( [4.121[ ) shows that 

which implies 

yt'^K* ^= ds<K*co^. (4.126) 



Since -ftT* = ^[1 + 0(1)], we obtain y^^^ < coet and thus f = oo. This shows 

e 



^ yf'^' ^ K*co-^ for A^^t^T. (4.127) 



3. Let T ^ ^/e and ^ yi < 2/2 ^ ^*(''") ~ ^(''") be given. Let and y^^^-* be solutions of 
( 4.I2ID with initial conditions y^^ = yi and y^^ = y2, respectively. Then there exists 
a G [0, 1] such that the difference zt = yf"^ — y[^^ satisfies 

= -dJix*it) - y^P - Oz, t) ^ -r?a(t)z, (4.128) 
where we have used ( [4.116| ) and (|1|). It follows that 

^ yf^ - yf' ^ {y2 - yi) e-''"(*'^)/^ (4.129) 

which proves ( [4. 118] ) in particular. If r ^ 4-y/e, we can use the relation x*(t) — x^^*'^ = 
yt^ + (yf - yf""^) to show that 

x*{t) - ^ ^*co^ + (x*(r) - x(t)) e-''"(*•")/^ (4.130) 

which proves (^ITtI) for r ^ 4^. Finally, if ^ ^ 4^, we can use the fact 

that x*{t) — ^ x*{t) — x^'^*'^^ to prove that ( 4.1171) holds for some constant 

C > 0. □ 
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Let us now consider the process yt = yj = xt—xf'^^''^, starting at time r in y,- = 0, which 
describes the deviation due to noise from the deterministic solution xf'^^'^ . It satisfies the 
SDE 

dyt = - [a^t)y + b^{yt, t)] dt + ^ dWt, (4.131) 



where we have introduced 

dot 



(4.132) 



a^{t) = dj{xr'\t) 
b^{y,t) = /(xf + y,t)- fixT'n - a^t)y. 

The following bounds are direct consequences of Taylor's formula and Proposition [4.11 



a*{t) ^ a^{t) ^ a{t) (4.133) 



a^{t) = a*{t) + + 0(te-^"(*'")/") (4.134) 

{a^Yit) = o(l + j e-^"(''^y^^ (4.135) 

|6^(y,t)| ^ 3My2(x*(t) + |y|), valid for x*{t) + \y\ ^ d. (4.136) 

For comparison, we will also consider the linear SDE 

dy° = -a^t)y^ dt + ^ dWf (4.137) 
e ye 

Let a'^{t, s) = a'^{u) du and denote by 

v^{t) = — f e^"" ds (4. 138) 
the variance of y^ . Again we introduce and investigate a function 

2|a(r)| e 7^ 

Lemma 4.12. T/ie function C^(i) satisfies the following relations for t ^ t ^ T : 



1 , . 1 



(C)'(t) ^ i. (4.142) 
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Proof: 

1. By integration by parts, we find 



at) 



1 



2\a{t)\ 



ds. 



(4.143) 



The relation |a'^(s)| ^ \a{s)\ ^ ?7|o(s)| together with (4.135) yields 



^ const + ^ e-'^"("'^)/^) e-2'?"(*'^)/^ ds. (4.144) 



The second term in brackets gives a contribution of order j e~v"(^''^y^ . In order to 
estimate the contribution of the first term, we perform the change of variables u = 
rj^t"^ — s^)/2e, thereby obtaining 



Vt' Jo (1 - n/e)3/2 ^ 7?t3 



23/^ + 2^!^ 



(4.145) 



where ^ = r]t'^/2e and = f/'7"^/2e. The last inequality is obtained by splitting the 
integral at ^/2. Using the fact that t^e'"'*^/'^^ < (6e/r/)3/2 ^-3/2 ^j. ^ ^ 

we reach 

the conclusion that this integral is bounded by a constant times e/t"^, which completes 
the proof of ( pIo| ). 
2. We now use the fact that C^{t) solves the ODE 



(4.146) 



Then, ( [4.142| ) is an immediate consequence of this relation, and ( [4.141| ) is obtained 
from the fact that 



dt 

whenever C^(^) = l/2|fi(''")|) and 



dtif' 2|a'^(t)| 



1 



KM 

\a{r)\ 

Jam 
\a^it)\ 



+ 1)^0, 



+ 1 



a*'{t) 
2a*(t)2 



> 0, 



(4.147) 



(4.148) 



whenever C^{t) = l/2|a*(r)|. Here we used ( [4.133[ ) and the monotonicity of d{t) for 

small t. □ 

We note that Lemma 4.12| and the bounds ( 4.133| ) on imply the existence of constants 
c+ ^ c_ > 0, depending only on / and T, such that 



- ^ at) ^ ^ 
t ^ ' t 



VtG[r,T]. (4.149) 

We can now easily prove that remains in a strip of width hy/^ with high probability, 
in much the same way as in Proposition 3.3. 
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Proposition 4.13. For sufficiently small T and e, and all t € [t, T], 

.01 . . ^^2 



Ws 



sup 



e) exp 



where r(e) = 0{e) and 



2^2 



+ 2. 



(4.150) 



(4.151) 



Proof: Let K = \\a^{t, r)|/2e2] and define a partition r = < • • • < = t of [r, i] by 



A; = 1,...,K-1. 



(4.152) 



Since a^(s) ^ a(s) ^ — r/s/2, we obtain — ^ Ae^ / {rjUk^i) for all A;. Now we can 
proceed as in the proof of Proposition |3.3|. □ 



We can now compare the solutions of the linear and the nonlinear equation. To do so, 
we define the events 



Vlt{h) = [uj: \yl\ < hy^Ci^ys G [r,t]} 



(4.153) 
(4.154) 



The following proposition shows that yl and differ only slightly. 

Proposition 4.14. Let 7 = 1 V 48M(2 + ^Jc^^c?^^ j ^fcZ and assume h < r/7 as well as 
hi^[d- x''{t)]^/T/{2^). Then 



n'^tih) c' fit 



1 + 7- 



1 + 7- 



h 

h]. 



Proof: Assume first that lo G Qf{h). We introduce the difference Zs 
6 = jh/r < 1, and define the first exit time 



(4.155) 
(4.156) 

Vs, set 



f = inf{s G [T,t]: \zs\ > 5hy^C{s) } G [r,t] U {00}. (4.157) 
On A = n^{h) n {f < 00}, we get by the estimate ( [4.136| ) on 6^, Lemma g?T| and ( |4.149D 



^ 6M(1 + (^)2 (2c+- + (1 + <5)4/'^) ^ /i7C(^ < <5/iycR, 



(4.158) 



for all s G [r, f], which leads to a contradiction for s = f. We conclude that ¥(A) = and 
thus \zs\ ^ jh'^ \/C'^{s)/t for all s in [t, f], which proves ( [4.156| ). The inclusion ( f4.155 ) is 
a straightforward consequence of the same estimates. □ 

Now, the following corollary is a direct consequence of the two preceding propositions. 
Corollary 4.15. There exists Hq such that if h < h^r, then 



T,t,x{t) 



}, (4.159) 



where C'^{t,e) is given by (4.151). 
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Appendix 

The appendix provides two lemmas needed in Sections ^ and ^. The first one uses expo- 
nential martingales to deduce an exponential bound on the probability that a stochastic 
integral exceeds a given value. 

Lemma A.l. Let (p{u) be a Borel-measurahle deterministic function such that 

ft 



<^>{t) = I ifiuYdu (A.l) 

exists. Then 







Proof: Let P denote the left-hand side of ( |A.2|) . For any 7 > 0, we have 

P = P| sup exp|7 / ip{u) dWu} ^ e^^l ^ pj sup M, ^ e^''-^*^), (A.3) 



where 



Ms = exp{^ jifiu) dWu-^J^ l\{uf dn} (A.4) 

is an (exponential) martingale, satisfying E{Mf} = E{Afo} = 1, which implies by Doob's 
submartingale inequality, that 

P| sup Ms-^ \] ^\^{Mt] = \. (A.5) 

This gives us 

P^e-T"^+^*W, (A.6) 

and we obtain the result by optimizing (A.6) over 7. □ 

The following lemma allows to estimate expectation values by integration by parts. 

Lemma A. 2. Let t ^ sq be a random variable satisfying Ft-{s) = P{t < s} ^ G{s) for 
some continuously dijjerentiable function G. Then 

W.{^s,,t){r)g{r)]^g{t)[F,{t)-G{t)\ + f g{s)G'{s)ds (A7) 

J so 

holds for all t > sq and all functions ^ g ^ 1 satisfying the two conditions 

• there exists an si G (sq, 00] such that g is continuously differentiable and increasing 
on (so,si); 

• g{s) = 1 for all s ^ si. 
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Proof: First note that for all t ^ si, 

ft , rtAT 



Pt f't/\T 

/ g'{s)F{T^s}ds = E{ g'{s)ds] 

J So J So 

= E{g{tAT)}-g{so) 

= E{l[,,,t)(r)5(r)} + g{t)F{T ^ t] - g{so) (A.8) 
which implies, by integration by parts, 

n\so,t){r)g{r)} = [ g'{s)[l-Fr{s)]ds-g{t)[l-Fr{t)]+g{so) 

J so 

^ f g{s)G'{s) ds + g{t) [Fr{t) - G{t)\ , (A.9) 

J so 

where we have used Fr{s) ^ G{s) and G{sq) ^ F{sq) = 0. This proves the assertion in 
the case t ^ si. In the case t > si, we have 

m{so,t){r)g{T)] = E{l[,„,,,)(r)5(r)} + P{r G [si,t)] 



r g{s)G'{s) ds + g{si) [Fr{si) - + [F^{t) - 

J so 

= f g{s)G'{s) ds - [G{t) - G{si)\ + [F^{t) - , (A.IO) 

J so 

where we have used that g{s) = 1 holds for all s G This proves the assertion for 

t>si. □ 
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